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Resolving Curvature Singularities in Holomorphic Gravity 
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We formulate holomorphic theory of gravity and study how the holomorphy symmetry alters the 
two most important singular solutions of general relativity: black holes and cosmology. We show 
that typical observers (freely) falling into a holomorphic black hole do not encounter a curvature 
singularity. Likewise, typical observers do not experience Big Bang singularity. Unlike Hermitian 
gravity 1], Holomorphic gravity does not respect the reciprocity symmetry and thus it is mainly a 
toy model for a gravity theory formulated on complex space-times. Yet it is a model that deserves 
a closer investigation since in many aspects it resembles Hermitian gravity and yet calculations are 
simpler. We have indications that holomorphic gravity reduces to the laws of general relativity 
correctly at large distance scales. 



I. INTRODUCTION 

The existence of singularities in Einstein's well tested the- 
ory of general relativity have puzzled many physicists 
since their discovery. It is a goal of (and motivation for) 
theories as quantum gravity to remove them. We propose 
generalizations of general relativity which to a great ex- 
tent ease the singularity problems already at the classical 
level of the theory. Singularities of general relativity are 
typically manifest as divergences of some curvature in- 
variant and are ubiquitous in general relativity ^, ^ . 

Consider, for example, the Big Bang singularity. A 
universe undergoing a power-law expansion expands with 
a scale factor, a ~ t^^"^, where e — —H/H^ ~ constant > 
denotes the 'slow roll' parameter, H — d/a is the Hub- 
ble expansion rate and dot denotes a derivative with re- 
spect to physical time t. The Ricci scalar R then diverges 
as t — + (corresponding to the time when matter density 
diverges) , 



6(2 -e) 

£2^2 



as t ^ 0. (1) 



Analogously, the Schwarzschild metric is singular when 
the coordinate radius r = ||a;|| — > (corresponding to the 
place where all of the mass is concentrated), resulting in 
the curvature singularity of the Riemann tensor, 



oo , 



as r 



0, (2) 



where M denotes the black hole mass, Gn the Newton 
constant and c the speed of light. 

On complex manifolds space-time coordinates x'^ get 
complexified as Q, ^ 



(3) 
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^ For the sake of simplicity here we drop a factor 1 / \/2 which we 
used in Ref. [J] when relating z'^ to and p^. 



where denotes the energy-momentum four-vector of 
noninertial frames, which can take any value; for an on- 
shell particle/observer however, pf^ reduces to the parti- 
cle's energy-momentum. Based on the holomorphy sym- 
metry, we expect that the singularities llHll) appear as 
some power of 
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(4) 



(and possibly its complex conjugate), where z ^ z'^ = 
J.Q _|_ jyO — _|_ i{G M / c^)E for the Big Bang singular- 
ity C]) and z ^ ||z || = ||x|| +z||y || = \\x\\ + i{G n / c^)\\p\\ 
for the black hole singularity ([2]). Near singularities the 
observed space-time curvature is typically proportional 
to a power of the real part of ((4|) , 



which blows up only when both x and y are simultane- 
ously zero. The aim is now to formulate a generalization 
of general relativity that yields such complexified singu- 
larities and show that for a freely falling observer x and 
y are almost never zero simultaneously. ^ In Hermitian 
gravity [l| , the Big Bang singularity can be considered as 
'resolved' in the sense that the set of observers moving 
backwards in time that encounter the Big Bang singular- 
ity is of measure zero. So typical observers do not see the 
singularity. Note that this is not the case in general rel- 
ativity, where all backward-moving observers eventually 
hit the Bing Bang singularity. In this paper we present 
a new theory, holomorphic gravity, in which also typical 
observers falling into a Schwarzschild black hole do not 
encounter a singularity. We believe that this resolution 
of singularities is a generic feature of gravity theories for- 
mulated on complex spaces as presented in this work and 
in Ref. [ij, representing one of the principal advantages 
of complex theories of gravity when compared with Ein- 
stein's general relativity. 



The meaning of 'almost never' is made more precise below. 
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II. ALMOST COMPLEX STRUCTURE 

A natural generalization of general relativity, in order to 
obtain complex solutions as (H]), would be to consider a 
theory with complex metrics, living on complex mani- 
folds. A general complex metric on a complex manifold 
is given by 



paper we construct a theory of gravity, based on the holo- 
morphic metric In the flat space limit the holomor- 
phic metric is invariant under the complexified Lorentz 
group, 50(1, 3; C), while the line element is invariant un- 
der the complexified inhomogcneous Lorentz group (the 
complexified Poincare group) /S'0(1,3;C). 



+Cu,,dz'' IS) dz" + Cniydzf" <S)dz^ 



(5) 



where barred indices = z^ denote complex conjuga- 
tion. Hermitian gravity [l| is formulated on a Hermitian 
manifold endowed with a Hermitian metric, defined as 
follows 



Cp {JpZ, JpW) = Cp (Z, W) , 



(6) 



where the action of the almost complex structure oper- 
ator, J, on the basis vectors of the complexified tangent 
space is given by 



''\^z^^)~ 'dz^- 



(7) 



A Hermitian metric is a complex metric which has - as a 
consequence of the symmetry requirement © - vanishing 
Cp,i, and Cpp components: 



C = C^pdz^ ( 



. dz'' Cn^dzf" ® dz" . 



The requirement that the complex metric satisfies 
Bianchi identities lead us in ^] to the conclusion that 
Eq. ^ can be consistently imposed on the metric only 
at the level of the equations of motion (on-shell) , while at 
the level of the action (off-shell) all complex metrics of the 
form ([5]) are in fact allowed. Thus studying holomorphic 
gravity allows one also to rectify the differences between 
the full complex theory (which, apart from holomorphy 
on vielbeins, has no additional symmetry requirements) 
and hermitian gravity. 

One of the reasons why one wants to study a theory 
of gravity that is invariant under the operation of J, is 
that the commutation relations of quantum mechanics 
are invariant under its action, which seems an invitation 
for quantization (recall that the coordinate in ^ is 
identified with the energy- momentum coordinate, p'', in 
the commutation relations). 

We can also consider the theory which is anti- 
symmetric under the action of the almost complex struc- 
ture operator, in the sense that the holomorphic metric 
is defined in the following manner 

Cp {JpZ, JpW) = -Cp (Z, W) . (8) 

The holomorphic metric can then be written as 

C ^C^.^dz'' ^dz" + Cf,^dzf ^dz"", (9) 

where the component C^^ (O^p ) is (anti-)holomorphic, 
which simplifies calculations. In the remainder of this 



III. THE HOLOMORPHIC METRIC 

The holomorphic line element is defined in the following 
manner 

ds^ = dz'^Cf.^dz" + dz'^Cf.i.dz'', (10) 
which can be written in eight dimensional notation 
ds^ = dz^ -C -dz^ {dz"^fC,nndz'' 

\ f dz'' 
dz"" 



(dz'^^dzf) 







(11) 



where the Latin indices can take the values 0, 1, — 
1,0, 1, c? — 1, the Greek indices run in the range 
0, 1, d — 1 and d denotes the complex dimension of the 
complex manifold. The entries of the metric C^n are 
functions of holomorphic and antiholomorphic vielbeins 
defined as follows ^ 



e{z)J'riabe{z)J' 



(12) 



Cf,u = e{z)^°-rjabe[z)^ . 

The metric components are symmetric under transposi- 
tion, since they are just an inner product of the viel- 
bein times its transpose. Hence also in eight dimen- 
sional notation the metric is symmetric under transposi- 
tion C'^ = C . We define the z^^ and z^ coordinates in 
terms of x^^ and j/'^, such that we obtain * 



d i / c .a 



1/ d 



d 



and their complex conjugates. This implies the following 
decomposition of complex vielbeins in their real, ej^J^, 
and imaginary, e/?, parts in the following manner 



< = - ieil ' ^ = = ej^^ + iei^. (13) 

Vielbeins are holomorphic functions, and thus transform 
as holomorphic vectors (we consider only the transforma- 
tion of the Greek indices for this purpose). 



dz'^ivj") 

el{z-)^~el{w-)-^^el{z^) 
e(^(z'')-eT(z.'^) = ^^^e?(z^). 



(14) 



^ Here the Latin indices a,b run from 0,1,..., d — 1, since they 
represent local indices, and r;^^ = diag(— 1, 1, 1, 1). 

* Checks ( " ) are put on indices to denote the imaginary part of 
a coordinate or on indices of objects, which are projected onto 
their basis vectors. 
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The holomorphy of vielbeins, 



IV. FLAT SPACE 



d 



d 



d 



\ a I - ai 



implies the Cauchy-Riemann equations, 



dyf' 



(15) 



When rotating the metric from z and z space to x and y 
space, the components of the holomorphic rotated metric 
are given by 



1 



(16) 



Its inverse then becomes 

1 / C^'^ + C^"' 



Clearly the entries of these rotated metrics are all sym- 
metric and real. It is easily verified that C^v = 9iiv—iKf^u 
and that Cpp = g^i. + iK^^, = 5^£> = gf,^ is a sym- 
metric metric tensor and gjxc = ~g^v We can write the 
metric in terms of the real and imaginary parts of the 
vielbein and in terms of the imaginary part of C^i/, using 
the definition of the complex metric in terms of vielbeins 
p2|. yields 



(17) 



Expressing the inverse metric gr™" in terms of vielbeins 
yields 



(18) 



With the rotated metric we can write down the holomor- 
phic line element in its rotated form 



ds^ = g^^dx^'dx" - ^^Kf^^dp^'dx" - ^gf.^^dp'^dp'' . 



From the definition of the complex metric in terms of 
vielbeins (I12|l and the fact that the vielbeins are defined 
to satisfy 

it follows that g"^"^ and g^n are inverses of each other 

g gen = o„ , 

where 5™ is defined as 

' Si^ 



(5™ = 



<5g 



The holomorphic line element p9)) in flat space becomes 

ds^ = -{cdtf + {dxf - [-{dyy + [dyf] . (20) 

Recall that the y coordinate is related to the energy- 
momentum coordinate [H, Q as = p^^Gn/c^- The 
space-time-momentum-energy interval squared from the 
origin to a space-time-momentum-energy point a;™ is 
given by 



d' {Q-x"') = ~{ctY + {x) 



G 



(21) 



Note that the contribution of momentum-energy to the 
line element multiplies G'j^/c^, which is tiny {Gj^/c^ ~ 
lO^^^s^/kg^), as it should be, since one does not observe 
any momentum-energy contribution to (|2ip at low en- 
ergies. For light-like propagation (ds^ = 0) and in the 
absence of momentum-energy contribution, Eq. (|2ip re- 
duces to the well known result: massless particles move 
on the light-cone with the speed of light, v = \\dx/dt\\ = 
c. 

On the other hand, for a light-like propagation in the 
presence of a non- vanishing momentum-energy contribu- 
tions however, Eq. ((2T|) yields. 



P^G 



N 



0, 



(22) 



where the variable p — 'q,_ivP'^p'^ is the energy- 
momentum squared in a space-time-momentum-energy 
hyper-surface. Setting the space-time-momentum-energy 
interval to zero determines the boundary of causality 
(light cones). We can specify our hyper-surface further 
by setting t to zero. Assuming p^ > 0, one finds that 
there is a spatial region which is in instantaneous causal 
contact, and whose radius is given by. 



/p^Gn 



(p2 > 0) . 



The causally related regions (time-like, ds^ < 0) for the 
(19) 

case p^ > are shown as white in figure (TJ where the 



region of non-local causal contact at t = is clearly seen 
as the hyperbola's throat. On the other hand, when 
p^ < 0, the causal boundaries, shown in figure [21 are 



quite different. There is a minimal time tmm 
for events to be in causal contact, 



required 



i 171. in, — 



N 



[p^ < 0) , 



We now calculate the phase velocity, using the space- 
time-momentum-energy line element 



^phase 



- _ - W ^2 



P^G^ 



N 



(23) 
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FIG. 1: The figure shows how the hght cones get modified 
by non inertial coordinate transformations on a space-time- 
momentum-energy diagram in the case when > 0. Note 
that there is a non-local instantaneous causally related volume 
element near the origin. 
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FIG. 2: The figure shows how the light cones get modified by 
non-inertial coordinate transformations when < 0. There 
is a minimum time interval required for events to be in causal 
contact. 

One can see that the phase velocity approaches the speed 
of light for large t. 

On the other hand, the group velocity becomes 

^9^\lc' + ^P, (24) 

where 




The group velocity approaches the speed of hght for small 
four-forces squared. Requiring that Vg be real and pos- 
itive (as required by propagation) sets a lower hmit on 
the four force, > —c^/G%. This is to be contrasted 
with an upper limit occurring in Hermitian gravity 
Just like in the case of Hermitian gravity, we expect that 
these aparent violations of causality become of impor- 
tance only in the regions where the particle's self-gravity 
is large (in the vicinity of particle's own event horizon). 

V. SECOND ORDER FORMALISM 

Let us now state the holomorphic Einstein-Hilbert ac- 
tion, ^ 

4 f 

S = / d'^zV^R + c.c. (25) 

IdttGat J 

The holomorphic Ricci scalar is defined by i? = C^'^R^^, 
where the holomorphic Ricci tensor is given by 

RUu = dx^l^. - 5,r^^ + r^r^^ - r^^r^^ , (26) 

and where here d\ = djdz^. Since the holomorphic 
Einstein-Hilbert action is just the complexified equiva- 
lent of the Einstein-Hilbert action of standard general 
relativity, the equations of motions for this theory are 
just complexified equivalents of the equations of motions 
of general relativity. Using first order formalism the 
equations of motion implied by the action (|25p are: 

Gpi/ = — — T^^l, (27a) 

VpC^'' = 0. (27b) 

Through the holomorphic metric compatibility equations 
()27bp the holomorphic connection coefficients are given 

by 

= \c'\dt.C,, 9.C^e - 5,C^,) . (28) 

Plugging the holomorphic Levi-Civita connection coeffi- 
cients into the holomorphic Einstein's equations (|27ap . 
we obtain second order differential equations in terms of 
the holomorphic metric only. 

When varying the holomorphic line element we 
obtain the following holomorphic geodesic equations 

F + V^^zn- = (29) 



^ In contrast to Hermitian gravity [H, the holomorphy symmetry 
of the metric tensor of holomorphic gravity can be imposed both 
on- and off-shell (at the level of the action). This has the advan- 
tage that the action II25I I suffices to fully specify the dynamics 
of holomorphic gravity, and no further (action) constraints are 
needed. Furthermore, it can be easily shown that the Bianchi 
identities and the covariant stress-energy conservation are satis- 
fied. 
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and its complex conjugate. We can write equations p9|) 
in their convenient eight dimensional form, such that we 
can rotate them to x, y space, by simply plugging in the 
rotated metric (|18p . This is possible because the holo- 
morphic connection coefficients transform as a (1,2) ten- 
sor under the constant coordinate transformations such 
as rotations from z, z io x,y space. The rotated eight di- 
mensional connection coefficients are then simply given 
by 



(30) 



where the metric g„j„ is just the eight dimensional ro- 
tated metric (fTBI). 



without any first order corrections present. We have not 
studied in detail the reduction of holomorphic Einstein's 
equations to general relativity. Our study of holomophic 
Schwarzschild solution and cosmology in sections IVIII 
and IVIIII suggests however that holomorphic gravity re- 
duces to Einstein's theory in the low energy limit. 



A. Scalar Field 

In this subsection we consider the following scalar field 
action 



1 



c.c.(,33) 



VI. THE LIMIT TO GENERAL RELATIVITY 

The limit of Holomorphic gravity to general relativity is 
based on the assumption that the y coordinate and its 
corresponding vielbein are small. When expanding the 
equations of holomorphic gravity in powers of y and its 
corresponding vielbein, we hope to obtain the theory of 
general relativity at zeroth order of the expansion and 
meaningful corrections at linear order. The easiest way 
to obtain the limit to general relativity is to expand the 
vielbeins in terms of the y coordinate 



el.{x,y) 



(31) 



From this and the definition, Cp,y = ^^(z) • e^{z), we 
conclude that the analogous formula holds for the holo- 
morphic metric tensor, C^i,{x,y) = e^^'^^C^vix), such 
that, 

g^uix,y) = coa{y ■ dx)g^^{x) +sm{y ■ dx)Kf_,i,{x) 
Kf^„{x,y) = - sin {y ■ da:) gf^^{x) + cos {y ■ dx) Kf_,^{x) , 

(32) 

where we made use of Eqs. ([TS]) and (|18p . The analogous 
relations hold for the inverses g^" and K'^'^. Assuming 
that K^i, is ~ just like - a first order quantity, Eq. ([5^ 
then tells us that g^^uix, y) acquires corrections at second 
order in y, while iir^^(x, y) acquires first order corrections 
in y, with K^^ being on its own a first order quantity. 

This analysis then implies that the connection with all 
indices unchecked yields the ordinary Levi-Civita connec- 
tion plus corrections of second order 

TP^Ax.y)^T%{x) + 0{y^). 

With these connection coefficients, we can know check 
if the theory reduces to the theory of general relativ- 
ity by plugging them into the rotated eight dimensional 
geodesic equation. Keeping only terms of linear order in 
the y coordinate and its corresponding vielbein, yields 
the ordinary geodesic equation 

xP + TP{x)x^'x'' + 0{y'') ^ Q 



where (j) = 4>{zP') is a holomorphic function of z^. The 
equations of motion for the scalar field are then 



(34) 



and its complex conjugate, where □ ~ C^V^V^ and 
□ — C'^'^V^Vp. When varying this action with respect 
to the metric we obtain the two sets of components of 
the stress-energy tensor 

= V„0V;50 - ^Ca^C'^'-V^^V.^ - Cc.fjV{(b) (35) 

and its complex conjugate. The energy conservation 
equations are 



and its complex conjugate, where we have used Eq. 
There are two sets of Einstein's equations, namely 



(V^0V,0 + C^,V{4>)) (36) 



and again its complex conjugate. The energy-momentum 
tensor of a perfect fiuid in the fiuid rest frame can be 
written in the diagonal form 



T^^ = diag(-p,p,p,p) 



(37) 



plus its complex conjugate, where p is the density and p 
the pressure. For a homogeneous scalar field = 0(t), 
Eq. dSni) yields, p = {Il2c^)4>'^+V{4)) andp = {1/2^)^'^- 

v{4>). 



VII. THE HOLOMORPHIC SCHWARZSCHILD 
SOLUTION 

The holomorphic Schwarzschild solution should just be 
a complexified version of the ordinary Schwarzschild so- 
lution, since the holomorphic Einstein's equations (|27ap 
are just complexified versions of the ordinary Einstein's 
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equations. So we expect to have the holomorphic metric 
components 



00 



Cii = 1 - 



1 2GnM 
c z 
1 2GnM 



(38) 



and their complex conjugates, where z ~ \\z\\ = 



i^^p and M is the mass of the black 



hole, which in holomorphic gravity could, in principle, 
be complex. ^ The angular part takes the same form as 
in general relativity 



dfl^ = de^ + sin' 



(39) 



but now the spherical symmetry is 0(3; C) and thus the 
angular parameters are complex 



= Or + tOi 



= 4>R + i(Pi- 



Thus the angular part of the holomorphic Schwarzschild 
solutions is the complex angular part ([55]) plus its com- 
plex conjugate. One can easily check that the holomor- 
phic Schwarzschild solution is indeed a solution to the 
holomorphic Einstein's equations ([27a|) in vacuiim. Ex- 
pressing the Schwarzschild components in r and p 
yields 



Coo — ^ 



2GnM r-i^p 



-I- 

2GnM 



(r + i^p) 



(r 



2GnM\2 



-7^(40) 

'-'OO 



and their complex conjugates. The rotated holomorphic 
metric components (jl6p are given in terms of the holo- 
morphic metric components (|40p . and take the following 
form 



900 = - 1 - 



2GnM 



Gi 



911 



9o6 



9ii 



„2 2GnKI 
■P TS r 



2GiM 



-P 



(r 



At the moment we do not have a good physical interpretation 
for S[M]. 



which reduces to the Schwarzschild solution of general 
relativity when p ^ 0. When these rotated compo- 
nents are inserted into the rotated holomorphic line el- 
ement ([T^ . one obtains the holomorphic Schwarzschild 
metric in its full glory 



2GnM 



Gl 



2GnM„ 



[cAtf 



^idEf 



(r 



2GnM 



2Gi 



{drf 



-P 



2G 



N 



(r 



2GjvM 



2Gi,M 



2Gn 



drdp 



A^rpdOndei 



dtdE 



f<2 



- dOA 



{d^R^ - dcbl^) X 



fi2 



p^ {{sin0RCOshdi)^ - (cos6iflSinh6l/)2} 



G 

— 4 — 5-rp(sin6'i^cosh(?/cos07jsinh6'/) 



2 

C° 



2sin6'fl,cosh0/cos0i^sinh6'/ 



- 2d(t)Rd(f)i 
G 

+ 2^rp{{smeRC0sheif - (cos6'7?,sinh6l7)2} 



This rotated holomorphic Schwarzschild metric is indeed 
a solution of the rotated holomorphic Einstein's equa- 
tions in vacuiim, as can be checked by explicit calcula- 
tion of Rfi^. It is easy to see that, in the limit when the 
radius r goes to infinity and p to zero, the solution ap- 
proaches the Minkowski metric, whereas the solution ap- 
proaches the momentum-energy Minkowski space when 
r goes to zero and p goes to infinity. If p is not zero 
when r is zero, there is no curvature singularity at the 
origin. Explicit calculation shows that p is not zero when 
r is zero for a generic infalling observer. When consid- 
ering an observer which is falling in radially we can ne- 
glect the change in angular coordinates, corresponding 
to a vanishing angular momentum. (Recall that in gen- 
eral relativity this choice corresponds to the worst case 
scenario, for which an observer descends in the quickest 
possible way towards the black hole singularity.) Based 
on the holomorphic Schwarzschild metric ((38)l the black 
hole line element pH)) of a radially infalling observer can 
be recast as 



~2c' 



1 



1 2GnM 



dz^ 
dr 



1 - 



1 2GnM 



dz \ 
J?) 



-r +C.C, 



(41) 



where we introduced a real affine parameter (proper 
time) T defined as c^{dT)'^ = —ds^. The holomorphic 
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Schwarzschild solution has an isometry in the direction 
{d/dz^ is a Killing vector) and hence there is a conserved 
quantity: 



-Itn + i- 



.G 



N fO 
JO 



(42) 



The real part = u"(to) corresponds to the (initial) 
zeroth component of the four- velocity evaluated at r = tq 
(the energy per unit mass divided by c) , which is also the 
Killing vector in general relativity, while the imaginary 
part /g = /°(ro) = (dp°/dr)(To) is the zeroth component 
of the observer's four-force. Since /q does not contribute 
in general relativity, the general relativistic limit should 
be obtained by setting /q — > 0, which is indeed the case. 
In analogy with general relativity, the conserved quantity 
can be written as Q 



dz^ 
dr 



1 2GnM 



(43) 



Equations (|4T|) and ([431) are not enough to fully deter- 
mine z = z{t), because nothing is known about the imag- 
inary part of the expressions appearing in Eq. (|4ip . Since 
we are in holomorphic gravity, it is reasonable to demand 
z = z(t) to be a holomorphic function. With this ana- 
lytic extension we can now completely determine z ~ z{t) 
by inserting the conserved integral ([^5]) into ([1T|) . to ob- 
tain 



dz 
cdr 



z 



J, z 



2GnM 



(44) 



and the complex conjugate of this equation. Note that 
is just the Schwarzschild radius and the on-shell value 
oi b = br + ibi can be expressed in terms of the initial 
3-velocity, &^ = 7^ - 1 = {'jovo/cf, -fo^ = 1 - {vo/cf. 
Solving l|44p for proper time r yields 



bi 



-In 



C 



From this it follows that, when b is complex and r ^ 0, 
p is not zero. Moreover, one can show that quite generi- 
cally, when r — > 0, p grows large, which limits the growth 
of the curvature invariant which for Holomorphic 
gravity has the simple generalization. 



^^l^LV per 



\2ri 



(45) 



0, this reduces to Ri^^paW^" 



For example, when r 
— 12c^^r^/(GArp)^, which is negative. 

In order to verify the limited curvature conjecture, we 
solve (|44p numerically. For initial conditions, where the 
imaginary part of b is nonvanishing, p is indeed nonzero 
when r = 0, as can be seen from figure 3(a) In this 




(a) 




(b) 



FIG. 3; The real and imaginary part and the absolute value of 
the complex radial coordinate 2: as a function of proper time 
T for a freely falling observer for different initial conditions 
h. 3(a) The real and imaginary part and the absolute value 
of z with a complex initial condition h. In this case |z| = 
is never reached, indicating that the observer encounters no 
black hole singularity. 3(b) The real and imaginary part and 



the absolute value of z for a real initial 6. In this case both '^\z\ 
and 9 [2] reach simultaneously zero for some value of proper 
time r, indicating a black hole singularity. 



case there is no curvature singularity. When 5 [6] = 
however, then r and p simultaneously go to zero, the 
geodesies end at a curvature singularity, just like in the 
case of general relativity. This situation is illustrated in 
figure 3(b) Even though numerical solution indicates 
that the evolution continues after r = = j> is reached. 



at that point there is a branch point of the evolution, and 
the numerical integrator picks one of the Riemann sheets 
(this is indicated by the cuspy feature of the numerical 
solution at r = p = in figure 3(b)). Moreover, at 
r — ^ — p the curvature invariant (|45p becomes singular, 
as can be seen in figure [S] Since 3[6] = only when the 
initial force /o = exactly, the set of initial conditions 
where b is real is of measure zero, when compared to the 
set of all initial conditions, where b can be an arbitrary 
complex number. Hence we conclude that the curvature 
singularity is not seen by most of infalling observers. 
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FIG. 4: The real and imaginary part and the absolute value 
of the curvature invariant (|45p . which we for simplicity denote 
in figures by . as a function of proper time t for different 
initial conditions h. An infinite curvature is attained only 
when = 0, which represents a set of initial conditions of 
Fig. ^ 



measure zero. 



The real and imaginary part and the 
absolute value of the curvature invariant for complex initial 
conditions h. The curvature invariant remains finite at all 



times. Fig. 4(b) A part of the upper panel 4(a) zoomed in. At 
large distances the real part of the curvature invariant {solid 
black line) approaches the general relativistic solution, but 
close to and within the Schwarzschild radius the curvature 
exhibits a very different behavior, becoming even negative. 
The deviation from the general relativistic behavior becomes 
large only within the Schwarzschild radius and hence remains 
hidden behind the black hole horizon. 



Finally, we need to interpret what it means for r to 
become negative, see figure [H In the view of Eq. , 
for instance, we can take the sign out of the denominator 
and put it in the numerator 



M 



-M 



-r -\- i{GN /c^)p r — i{GN/c^)p 



(46) 



In this way we effectively glue the holomorphic 
Schwarzschild solution with a negative ^[z] onto the anti- 
holomorphic solution with a positive ^[z] but with a neg- 
ative mass: a white hole. This means that an observer 



FIG. 5: The real and imaginary part and the absolute value 
of the curvature invariant (|45|) with initial conditions corre- 
sponding to a real b (vanishing initial force). Since Q[z] = 
remains zero, the curvature invariant is identical to the cor- 
responding curvature invariant in general relativity ((T} such 
that a curvature singularity is reached at r = = p. 



can fall through a black hole and emerge from a white 
hole with an opposite momentum and v. v. Needless to 
say, this behavior is very different from general relativ- 
ity, where classically nothing can come out once inside 
the event horizon. Finally, we note that the holomor- 
phic and anti-holomorphic solutions are in general not 
symmetric under time reversal. This symmetry is bro- 
ken, for example, by the imaginary part of the curvature 



invariant, as can be seen in figure 4(a) 



VIII. COSMOLOGY 

Since our Universe is isotropic and homogeneous on large 
scales, our complex generalized theory should also pos- 
sess isotropy and homogeneity such that it reduces to 
the theory of general relativity correctly at low energies. 
Since our complex theories are completely specified in 
terms of vielbeins, an assumption concerning the viel- 
bein modeling an isotropic and homogeneous universe is 
in place here. Let us try the following Ansatz: 



(47) 



where 



a(z°) = flflX^c) + *a/(^c) «(^c) = aR{z°) - iai[z°) , 

where z^ = rj -\- iG^Ec denotes a conformal complex 
'time.' (In this section we set the speed of light c = 1.) 
With these assumptions the holomorphic Christoffel sym- 
bols (|28p are then given by 



a 
a 



(48) 



and its complex conjugate, where a' = da/dz^. Inserting 
Eq. into the holomorphic Riemann tensor we 
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obtain 



+ 



a/^_ fa/_ 
a \ a 

2 

[2SX + 2ry,,. 



(7?^. - 2SX) 



and its complex conjugate. The holomorphic Ricci 
scalars are then given by 



i? = 6— = 6 



and its complex conjugate, where we have contracted 
holomorphic Ricci tensors with the inverse metric g'^'^ = 
(l/a^)?]'"' and its complex conjugate, and d — da/dz^, 
dz^ = adz^. The holomorphic Einstein's equations in 
four dimensions are 

R^^l' — SnGNiTfiu — -g^i/T) 

and its complex conjugate, of course. Because of the 
isotropy and homogeneity symmetries imposed, the Ein- 
stein's equations contain only four independent equations 
due to homogeneity and isotropy symmetries, namely one 
space-space (ii) equation. 



a V a ,' 



the time-time (00) equation, 
AttG 



a 
a 



{p + ip) 



(49) 



(50) 



and the complex conjugates of these equations. We can 
simplify Eq. ([^5]) by making use of the time component 
(1501) to obtain 



SttG 



N 



(51) 



where H — d/a denotes the holomorphic Hubble parame- 
ter. Equation ((5T|) together with equation ((50)) constitute 
the holomorphic Friedmann equations. It is easy to see 
that the holomorphic Friedmann equations are consistent 
with the holomorphic energy conservation equations 



p + 3-ip + p) 
a 







(52) 



and its complex conjugate, which is derived from the 
holomorphic energy conservation equations 

and its complex conjugate, where we have used the ex- 
pression of the energy-momentum tensor (|37p and holo- 
morphic Christoffel symbol ([35]). Making use of Eqs. ([57]) 



and (j35p . we can express the energy density p and pres- 
sure p in terms of the scalar field (j). In an isotropic and 
homogeneous universe the scalar field is a holomorphic 
fimction of z*', (j) = (piz'^), such that we have 



P = l<i^' + v{^) 



p^l^2^V{cl>). (53) 



The equations for p and p are just the complex conju- 
gates of these expressions. The equation of motion for 
the scalar field (IMl) becomes 



• dV 
3iJ0 + ^ = 



(54) 



and the holomorphic Friedmann equations (|50H5ip then 
become 



SnG 



N 



a 3 



and 



a \ SttGn ' 



3 \ 2 



(55) 



(56) 



and the corresponding complex conjugates. Just like in 
general relativity, one can show that only three equations 
among Eqs. (|5^ and ([5iH55|) are independent. 

We shall now study the holomorphic Friedmann equa- 
tions for different physical circumstances and compare 
the results with those of general relativity. 



A. Power law expansion 

Let us now consider a homogeneous and holomorphic 
fluid with the equation of state, 

= weC, ^[w]>-l. (57) 

P 

In this case Eq. ([5^ is solved by, 

Next Eqs. ([55]) and (|56p can be combined into, 
H d / 1\ 3, 

such that e e C is a complex constant and 5R[e] > 0. 

When equation (|59p is integrated, one gets the holo- 
morphic expansion rate. 



H 



(60) 



This is easily integrated, to yield the holomorphic scale 
factor, 



Co 



(61) 
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These solutions are a generalization of power law expan- 
sion of general relativity, in that e is in general a complex 
parameter, is an arbitrary constant parameter that 
signifies the complex time at which a = 1. Its physical 
significance is revealed by realizing that Co parametrizes 
the Hubble rate at the time z'^ = (q, as can be seen from 
Eqs. PI | and H ^ . 

Finally, the last equation we need to solve is 
which gives. 



1 



SttGnPo 



(62) 



where Hq = i?(Co), Po = p(Co) 



We shall now show that, non unlike in general rela- 
tivity, the holomorphic power law expansion (|6ip can 
be realised by a homogeneous holomorphic scalar field 
(f> = (piz^) with an exponential potential. Just like in 
general relativity, the power law solution is then realised 
in the scaling limit, in which case it exhibits attractor 
behavior 

Let us consider the following Lagrangian density 



M 



(63) 



We now assume that the Friedmann equation permits 
a scaling solution of the form, cx l/z°. Obviously, such 
a solution must satisfy. 



Co 



z" 



ooCo In ( -y- 



(64) 



where 4>o,(t)o and Co are (complex) constants (note that 
Co and 00 ai'e not independent; indeed, rescaling Co can 
be compensated by the appropriate shift in (jjo). Note 
that Eq. (pT|) implies that the scalar field can be 
considered as a 'marker' for the scale factor. 



a = exp 



£00 Co 



(65) 



and in this sense defines a clock. 

A constant w requires the scaling of the potential ((63)) . 



T/ T/ I ^'l^o\ f Co 



(A0oCo)/A/ 



(66) 



We are free to absorb 00 in Vo by redefining, 
Vo exp {—X(f)o/M) Vo- Demanding the scaling, 



V^V,[^^] , (67) 



implies the following condition on the initial field velocity, 

2M 



ACo ' 



(68) 



where yo=n<^o) = n'/'(Co))." 

Taking account of these scalings, the equation of state 
parameter w becomes indeed constant. 



^ - 1 



2M^ 



r;o + l' 2^0 X^V^Ql 

and hence e in (|59p is also constant, 
d ( 1 



^ 3??o 
dz^ yn J ^ 1 + rio 



(69) 



(70) 



The conservation equation is now trivially satisfied 
implying the energy density scaling, p — po/a^'^ = 
Po{C,ol z^Y ■ Finally, the Friedmann equation (j56[ 
and (|62p yields the constraint 



2 87rGArVo 

tlr. = 

° 3-e 



(71) 



When this is combined with Eq. ([70]) . one gets the fol- 
lowing algebraic equation for e 



e(6-3) 6 



A2 



167rG'ArM2 



= 0, 



(72) 



which links the parameters of the scalar theory with e. 
Note that e does not depend on the initial conditions 
on the field, but only on the coupling parameters of the 
potential. This should not surprise us since we are con- 
sidering an attractor solution. 

From Eq. ([72]) it follows that de Sitter space (e — 0) 
and kination (e = 3) are solutions and have a special 
significance. One can show that 5R[e] < 3, the limit is 
saturated when = \&ttGisiM'^. For larger values of 
A there is no scaling solution: V redshifts faster than 
the kinetic term, resulting in e — > 3. de Sitter space 
e = is also a solution of Eq. ((72)) . which is realised when 
00 = 0. This is indeed a stable solution, for which kinetic 
energy vanishes, and V = Vq = Ao/{8ttGn), where Ao 
denotes the equivalent cosmological constant, which in 
holomorphic gravity can be complex. 

The nontrivial solution of Eq. ([72|) is 



A2 



< 5R[e] < 3 



(73) 



which corresponds to the attractor solution when the con- 
dition < di[e] < 3 in Eq. ([73]) is satisfied. Of course, 
the solution (j73p is more general than the correspond- 
ing solution of general relativity in that A and hence also 
e can be complex. To appreciate the significance of a 
complex e, let us assume that the Universe follows the 



One can show that, even when the condition II68I I is not met, 
quite generally ij> eventually approaches the attractor solution for 
which ^ holds. 
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attractor behavior with e given by (|73p . In this case the 
exponential potential (|63p can be written as 



V ^ Vq exp Vl67rGAre 0j , < 5R[e] < 3 . 

(74) 

This potential can be used to obtain both an accelerating 
universe (for which 3fi[e] < 1) and a decelerating universe 
(with 3?[e] > 1). Thus with the appropriate choice of 
A all standard cases in cosmology can be reproduced: 
radiation era (e = 2); matter era (e = 3/2); kination 
(e = 3), which is realized in the limit when A — > and 
Vb -> 0; inflation (0 < 3fi[e] < 1), etc. 

The physical Hubble parameter TL is given in terms of 
the real part of the holomorphic expansion rate H as, 
rC^ = For a real e Eq. ^ imphes, 



j-2 _ p2 



(75) 



For a complex e the expression is more complicated. 



(76) 



In order to find out whether the holomorphic Big Bang 
singularity is ever reached, we need to study geodesic 
equations, which shall give us a crucial information on 
whether a freely falling observer experiences the holo- 
morphic singularity in ([75H76)) (c/. Ref. [Ij]). 

The physical scale factor A can be obtained as the real 
part of the holomorphic scale factor squared. 



^,_,t' + G%E^ 



^--^Arcta 



X cos 



ICoP 
26/ , 2eR 



Arctan 



X In 



GnE 

e + GiE^ 



ICol 



(77) 



where e = e/j+ie/. When e/ = and provided E does not 
grow with time (which is reasonable), the Universe ap- 
proaches the standard FLRW cosmology. When however 
£/ 7^ 0, the Universe's scale factor develops oscillations, 
which can result in significant differences between holo- 
morphic and FLRW cosmology even at late times. This 
is a disadvantage of holomorphic cosmology when com- 
pared with, for example, Hermitian cosmology developed 
in Ref. [l|. Note that, when e develops an imaginary part, 
then the potential ([71]) violates charge-parity symmetry. 
In this (/) evolves, 3?[y] oscillates and can be ei- 

ther positive or negative. When the potential is negative, 
the Universe can enter an anti-de Sitter-like phase. As 



Note that the scalar potential of hermitian gravity discussed in 
Ref. Q differs from Eq. | [74t . 



a consequence, the physical scale factor J? in Eq. ([77|) 
can be either positive or negative. To prevent a negative 
value for one can add a constant to V (or a cosmolog- 
ical term), which will keep A^ positive at all times. This 
type of behavior can have relevance for the Universe's 
dark energy. 



B. Geodesic equation 

In order to better understand the behavior of the ex- 
pansion rate and the corresponding scale factor, we shall 
now consider a freely falling observer in the contracting 
phase. To do that, we need to solve the correspond- 
ing geodesic equation, which in holomorphic gravity has 
formally the same form as the corresponding geodesic 
equation of general relativity discussed for example in 
Ref. T|. Taking account of the Christoffel symbol (|48p . 
the geodesic equation and the line element can be written 
in conformal time as. 



6^. 



0, 



(78) 



a^77Q/3M"uf -I- a^rjafiu'^u^ = -2 



where t is the (real) proper time of a freely falling ob- 
server (in the frame in which all 3- velocities vanish): 
(ds)2 = -2{dTf, and 



dT 



(79) 



is the 4-velocity in conformal coordinates = (77,2;*). 
Defining the physical 4-velocity as 



au': 



'17 



(80) 



we can rewrite the spatial and time component of the 
geodesic equation ([78]) as. 



d{au^) 



0, 



d{a^[{i 



1]} 



= 0. 



dr dr 
These are solved by the following scaling solution. 



0\2 



(uO) 



1 



oo 
a 



(81) 



(82) 



where 



'0 



it^(To) and m'' = u'^{t). We introduce a 



complex constant U by recasting the first equation in 

as 



± 



U 



df_ 



Ur + lUi 



(83) 



When the line element ((78|) is taken account of, one would 
be tempted to identify. 



U 



al\\uof . 



(84) 
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From the line element ([78]) it follows that, strictly speak- 
ing, only the real part of Eq. ([55]) must be satisfied. Yet 
by holomorphy (which we assume to dictate a unique an- 
alytic extension) we know that also the imaginary parts 
must match. We assume here that Eq. (|83)) indeed holds 
true. 

Upon rewriting Eq. (j83p in the integral form we get. 



daa'- 



T 

<0 



(85) 



where we made use of Eq. ((6T|) and dz^ = ez'^da/a = 
eC,Qa'^~^da. Integrating ([55)1 yields a hypergeometric 
function. Rather than performing a general analysis of 
Eq. ([55]) , we shall consider the cosmologically interesting 
cases which are simple to analyze. We shall first consider 
a curvature dominated epoch with e = 1 . 



1. Curvature Dominated Epoch 



In this case Eq. ([85|) integrates to. 



(86) 



Note that that is not a unique expression; there is a 
freedom to shift r for an arbitrary real constant. From 
Eq. JHS]) it follows that, 



(87) 



= 

from where we conclude, 
1 




1 



From ((87)) we get for the physical scale factor. 



A' 



Wo 



M oe 



ICoP 
, Co 



cos(20J - \u\lcos{2e^ 
^iCole^'^S 



(89) 



This then implies that the scale factor A'^ will be positive 
at all times, provided both cos(26'ij) > and cos(26'„) < 
0, in which case the Big Bang singularity will never be 
reached. In order to have a more precise understanding 
on whether the Big Bang singularity is ever attained, we 
need to study a curvature invariant, one example being 
the Hubble parameter. 

Observe that the physical Hubble parameter, which is 
obtained from Eq. 



"|glCopcos(z9) 



|^|2|CoPcos(,9))2 + |^|4|Co|4sin^(^) 



(90) 




(a) 




(b) 



FIG. 6: The physical Hubble parameter (|90|) in the curvature 
dominated epoch as a function of the proper time for different 
values of •& remains finite at all times. An infinite value of the 
physical Hubble parameter is attained when Sjpo] = 0, which 
represents a set of initial conditions of measure zero compared 
to the set of all initial conditions. 



6(a) The physical Hubble 
parameter with initial conditions cos (W) < 0. [6tb)] The phys- 
ical Hubble parameter with initial conditions cos(i9) > 0. 



represents a bounce universe, where we defined i9 — 
2{du + 0c)- When cos{d) > the maximal physical 
expansion rate is reached when = |it|g|Cop cos(i^) -I- 
|u|olCoP sin(i?), for which 



2|^|^|CoPsin(^) 



(91) 



As can be seen in figure 6(a) 

^2 



this corresponds to a 
local maximum. At = |w|oTCoP cos(??), 7i = 0; at 
even smaller proper times < 0, which means that a 
local observer will have the impression that the Universe 
has entered an anti-de Sitter-like phase. Since we are in 
holomorphic gravity, there is no need to change the form 
of the line element (vielbein) Ansatz (HT)) . The minimal 
expansion rate squared 7i^i„ is reached when t = 0, 



cos{d) 



(92) 
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which is singular only when uq = 0, or equivalently when 
Mq = 1, corresponding to a set of initial conditions of 
measure zero. 

It is interesting to note that when cos(i?) < 0, the 
expansion rate (|92p becomes a global maximum, away 
from each Ti. decreases monotonously in both directions, 
as can be seen in figure [6(b)| This case represents a more 
conventional bounce Universe, and it is realized when 
the initial 3-force dominates over the initial 3-velocity, 
G^||/(ro)ll>IKro)||. 



2. Radiation Era 

Let us now consider radiation era (e = 2), in which case 
Eq. ([55]) integrates to. 



a\/ d?- + [/— C/ln (a + \J a? ^U] = — (radiation era) . 
^ ^ Co 

(93) 

This transcendental equation cannot be solved for a = 
a^r) in terms of elementary functions and thus it is hard 
to analyze in complete generality. A rather conclusive 
analysis can be, nevertheless, performed by integrating 
Eq. ([83|) numerically, which shows that, when U and 
Co are chosen real, the Hubble expansion rate — 
l/(2z")^ = l/(2Coa^)^ becomes singular at r = 0, as can 
be clearly seen in figure 7(a) WhenU — (M°)n— 1 = llulln 



is complex however, then ~ as can be seen in figure 7(b) 



both 3fi[i/^] and are in general nonzero, implying 

that the physical Hubble parameter remains finite at all 
times. We conclude that, just like in the case of curva- 
ture domination, there is no curvature singularity except 
for a set of initial conditions of measure zero (5[uo] = 0), 
when compared to the set of all initial conditions (3[ito] 
arbitrary). The Big Bang singularity is hence resolved 
also in the radiation era. We think that an analogous 
conclusion holds for more general expanding space-times. 
Even though the physics of radiation era is quite different 
from that of Schwarzschild black holes, the correspond- 
ing geodesic equations (|44|l and ([83]) - based on which we 
performed the analyses of singularities - are of an iden- 
tical form, and thus the conclusions of the analyses are 
quite similar. 



IX. CONCLUSIONS 

We have shown that - just like Hermitian gravity [ij 
- holomorphic gravity also resolves the cosmological Big 
Bang singularity. 

Furthermore, we have shown that holomorphic gravity 
resolves the Schwarzschild space-time singularity. Given 
the similarities between the Hermitian and holomorphic 
gravity theories, we expect that Hermitian gravity also 
does not suffer from the singularity problems of general 
relativity 

Even though, strictly speaking, singularities in com- 
plex theories of gravity exist, they occupy just one point 
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FIG. 7: The physical Hubble parameter in radiation era 
(e = 2) as a fu nction of proper time r for different initial 
conditions. 7(a) The physical Hubble parameter for observers 
with real initial condition U defined in Eq. (|83|l 



becomes sin- 
gular at some time t, which corresponds to the Big Bing sin- 
gularity. 



7(b) When U in Eq. (|83p is complex, the physical 
Hubble parameter remains finite at all times, thus resolving 
the Big Bang singularity. The Hubble parameter is in general 
nonsymmetric under time reversal. Only when U is purely 
imaginary, the physical Hubble parameter is symmetric un- 
der time reversal. 



on an 8 dimensional phase space, and consequently they 
are reached only by special geodesies/observers. Indeed, 
we have argued that the space of initial conditions of ob- 
servers which encounter singularities is of measure zero, 
when compared to the space of all initial conditions. 
Hence quite generically observers will not see singular- 
ities within holomorphic gravity. We expect that quanti- 
zation of space-time and momentum-energy will lead to 
further smearing of these singularities. 
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